Running phonon dynamicsΒΆ
In this last step, the resulting atomic dynamics are computed, with time taken into consideration. The dynamics are computed under 60000 time-steps. Phonon-phonon evolution are modeled using the below equation. The first portion of the equation models absorption and the latter emission. Delta_t here is infinitesimal (0.0005e-12).
For electron-phonon and hole-phonon evolution,
the below equation is used. Again, delta_t here
is infinitesimal (0.0005e-12).
After each time-step, the occupation states each of the particles are recalculated. The phonon occupation state changes due to the cumulative interactions it has with holes and electrons. Each hole and electron changes occupation states based on its temperature. As both particles either lose or gain energy in the process, temperature change can be found based on both particles' specific heat.